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Finite-temperature phase structure of lattice QCD witli tlie Wilson quark action is analyzed. We show that 
the critical line at finite temperatures, defined to be the line of vanishing pion screening mass, turns back toward 
strong coupling, forming a cusp on the (/3, K) plane, and that the line of thermal transition runs past the tip of 
the cusp without touching the critical line. Previous results are discussed in the light of our findings. 



1. Introduction 

Despite extensive effort over the years [^]-^ the 
finite-temperature phase structure of lattice QCD 
with dynamical Wilson quarks has not been fully 
clarified to date. The difficulty, originating from 
an explicit breaking of chiral symmetry with the 
Wilson action, arises from the question in what 
sense the concept of the critical line of vanish- 
ing pion mass at zero temperature extends to the 
finite-temperature system and how the line be- 
haves in the space of the inverse gauge coupling 
constant (3 = 6/17^ and the hopping parameter 
K. This is an issue of fundamental importance, 
not answered by previous studies, which only ex- 
amined the relation between the line of thermal 
transition and the zero-temperature critical line. 
In this article we report our resolution of the is- 
sue based on the idea of spontaneous break- 
down of parity and flavor symmetry 

2. Phase diagram for two flavors 

We define the critical line Kc{(3) for the finite- 
temperature system on a lattice of a temporal 
lattice size Nt by the vanishing of pion screen- 
ing mass extracted from the pion propagator for 
a large spatial separation. The line thus de- 
fined represents a natural extension of the zero- 
temperature critical line in that (i) it smoothly 
converges to the latter as Nt — > cxo, and that 
(ii) it corresponds to a singularity of the finite- 
temperature lattice QCD partition function. Let 
us also define Kt {(3) to be the line of thermal tran- 
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Figure f. Phase diagram for Nf = 2 with the 
Wilson quark action on an A^t = 4 lattice. 

sition, the thermal line. With these definitions 
the phase diagram we find for the case of Nf — 2 
flavors on an A't = 4 lattice is presented in Fig. |^. 
The solid line is our estimate of the critical line 
obtained by an interpolation of points where the 
pion screening mass vanishes (solid circles) or the 
number of conjugate gradient iterations in hybrid 
Monte Carlo runs is estimated to diverge (open 
circles). The dotted line is the thermal line with 
open squares representing its location reported in 
previous simulations 

Let us emphasize the novel features of our 
phase diagram. Contrary to the zero-temperature 
case, the critical line does not extend toward ar- 
bitrarily weak coupling, but rather turns back to- 
ward strong coupling at a finite value of /3, form- 
ing a cusp. The location of the critical line before 
the turning point is generally shifted from that 
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Figure 2. Critical line for the Gross-Neveu model. 
Lines correspond to temporal lattice size Nt = 
2, 4, 8, 16 and oo from inside to outside. 

at zero temperature. For an A^f = 4 lattice the 
shift is small, however. We expect the tip of the 
cusp to be rounded if examined at a fine scale, 
and that there will be four more of such cusps in 
the entire (/?, K) plane. The region bounded by 
the critical line represents a phase in which parity 
and flavor symmetry is spontaneously broken [^. 

We point out that the absence of the critical 
line toward weak coupling provides a natural ex- 
planation of the fact, observed in previous simu- 
lations |^,||, that physical quantities in the high- 
temperature phase vary smoothly across the zero- 
temperature critical line. 

Another important feature of our phase dia- 
gram is that the thermal line runs past the tip 
of the cusp of the critical line and continues to- 
ward larger values of K. In other words the low- 
temperature phase goes around the cusp and ex- 
tends into the region bordered by the upper part 
of the critical line. This contrasts with the con- 
clusion of the QCDPAX Collaboration that the 
two lines meet at /3 = 3.9 — 4.oQ. This point is 
discussed later. 

3. Evidence for the phase diagram 

3.1. Gross-Neveu model in two dimensions 

An instructive hint on the nature of the criti- 
cal line is provided by the two-dimensional Gross- 
Neveu model formulated with the Wilson quark 
action In the large N limit the critical line 
of this model can be calculated from the saddle 
point equations for the pion mass. The result is 
shown in Fig. on the plane of the bare coupling 
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Figure 3. Hybrid Monte Carlo results for parity- 
flavor order parameter {ipi^^T'i'ip) for the Nf = 2 
case at /3 = on a 4^ lattice as compared to 
predictions of \/Nc expansion (drawn by lines). 

constant g and bare quark mass m = 1/2K — 2, 
where the temporal lattice size is taken to be 
Nt = 2, 4, 8, 16 and oo from inside to outside. 

At zero temperature there are multiple criti- 
cal lines forming three cusps at g = 0. Outside 
of the critical line is the normal phase, while in- 
side is the parity-broken phase characterized by 
("075?/') 7^ 0. Pion is interpreted as the zero mode 
of the second-order transition separating the two 
phases 1^]. The conventional continuum limit is 
taken near {g,m) = (0,0) along the upper crit- 
ical line. The existence of two additional cusps 
arises from an exchange of doublers and physical 
fermion as m is varied. 

At finite temperatures the cusps move away 
from the weak-coupling limit g = 0. As a result 
the critical line forms a single continuous line. 

3.2. Zero-temperature simulations 

The chiral properties of the Gross-Neveu model 
are similar to those expected for lattice QCD. 
This leads us to anticipate that the structure and 
characterization of the critical line are also simi- 
lar except that there will be five cusps for QCD 
due to an increase of dimensions from two to four. 

The existence of a parity broken phase in QCD 
was originally suggested by an analytical calcu- 
lation in the strong coupling limit in the 1/Nc 
expansion|^. In Fig. ^ we show a comparison of 
the analytical prediction and our hybrid Monte 
Carlo results for the parity-flavor order parame- 
ter (ipij^T^ip) for the Nf — 2 case in the presence 
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Figure 4. Zero-temperature tt and p mass at (3 — 
6.0 in quenched QCD. 



of a symmetry- breaking external field H. A good 
agreement between the two results supports the 
conclusion of the large Nc analysis that parity 
and flavor symmetry become spontaneously bro- 
ken for K >1/A at f3 = 0. 

In Fig. ^ we plot the zero-temperature tt and 
p masses at f3 — 6.0 in quenched QCD as a func- 
tion of 1/K. The right-most set of data points 
are taken from the APE CoUaborationQ, while 
the rest of data are obtained with 10 — 20 con- 
figurations on a 10'^ X 20 lattice. In addition to 
the conventional critical line at 1/K = 6.3662, 
we observe clear indictions for three more values 
of Kc- While the fourth value expected toward 
^ is yet to be confirmed, this provides ev- 
idence that there are multiple critical lines toward 
weak coupling in QCD. 

3.3. Finite-temperature simulations 

In order to determine the shape of the finite- 
temperature critical line, we have carried out hy- 
brid Monte Carlo simulations for the Nf — 2 sys- 
tem on an 8"^ x4 lattice in the region 3.0 < /3 < 5.3 
and 0.15 <K< 0.31. 

At /3 = 3.5 results for the pion screening mass 
demonstrate that there exists a second critical 
value at Kc — 0.2454 beyond the conventional 
value at K^ = 0.2267 (see Fig. 2(a) in Ref. @). 
The behavior of the pion mass is similar to that 
around 1/K w 6.3 in Fig. ^. A parallel analysis 
at /3 = 4.0 shows that the gap between the two 
critical values is either extremely narrow or dis- 
appears just below /3 = 4.0 (Fig. 2(b) in Ref. g). 
Since previous simulations [|l]-^ found no singu- 



lar behavior beyond the thermal line down to 
(3 — 4.3[||, we conclude that the critical hne turns 
back toward strong coupling in the vicinity of 
/3 « 4.0. 

We expect the number of conjugate gradient 
iterations ncG for evaluating quark propagators 
in hybrid Monte Carlo runs to diverge toward the 
critical line. Extrapolating 1/ncG linearly in (3 at 
a fixed K, we estimate the location of the critical 
line for several values of K. The results, together 
with those from the pion mass, are smoothly con- 
nected, forming a cusp as is shown in Fig. |^. 

A characteristic feature of the high tempera- 
ture phase is that the pion screening mass takes 
a finite value. Such a phase can not include the 
critical line where the screening mass vanishes. 
Thus the region close to the entire critical line 
has to be in the low temperature phase. The va- 
lidity of this argument is confirmed through mea- 
surements of the quark and gluon entropy and 
the Polyakov line. Results at /3 = 3.5 show that 
these quantities decrease toward the second criti- 
cal value when it is approached from larger values 
of K^ and become comparable in magnitude to 
those below the conventional critical value which 
is in the low temperature phase (see Fig. 3(a) in 
Ref. 1). 

This means that the thermal line can not cross 
the critical line. Since the two lines are known to 
come close around /S ~ 4.0||^, we conclude that 
the thermal line runs past the tip of the cusp and 
continues toward larger values of K. Our results 
for thermodynamic quantities taken around the 
tip of the cusp at /3 w 4.0 (Fig. 3(b) in Ref. [|) 
and those in the literature [^js] are consistent with 
this conclusion. 

4. Discussions 

4.1. Comparison with previous results 

A naive expectation for the chiral transition 
with dynamical Wilson quarks is that the thermal 
line would intersect the zero-temperature critical 
line, which would naturally define the point of 
chiral transition in that it would be the point 
of vanishing pion mass separating the low- and 
high-temperature phases. Early work||l[] guided 
by this idea failed to find the crossing point. More 
recently the QCDPAX Collaboration concluded 
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that the crossing point with the expected prop- 
erties is located at /3 « 3.9 — 4.0 for the Nf — 2 
system on an iVt = 4 lattice^. 

However, the actual structure differs from their 
conclusion in one crucial aspect: while the ther- 
mal line does cross the zero-temperature critical 
line since the latter continues toward the weak 
coupling limit, the pion screening mass does not 
vanish at this crossing point. In fact a unique 
point marking chiral transition as naively ex- 
pected does not exist in our phase diagram. The 
whole region around the tip of the cusp and the 
thermal line is relevant for understanding the chi- 
ral properties of the thermal transition. 

4.2. Dependence on quark action 

The key ingredients of our analyses are the 
properties of the Wilson quark action that mass- 
less free quarks appear at multiple values of 
the hopping parameter and that the action has 
parity-flavor symmetry. We therefore expect the 
qualitative structure of our phase digram to re- 
main valid for a more general form of Wilson-type 
actions including the clover action. 

We remark that our phase diagram does not 
explain the metastability signals at /3 « 5.0 found 
by the MILC Collaboration|3j. This phenomenon 
is probably a lattice artifact closely related to the 
specific form of the Wilson plaquette or quark 
action 

4.3. Dependence on Nf 

We expect our phase diagram to apply also for 
a larger number of flavors, at least up to Nf = 6, 
except for one difference. While simulation re- 
sults strongly indicate that the thermal transition 
for the Nf — 2 case is a continuous crossover and 
not a true phase transition, first-order signals are 
found for for the case of Nf = 3 and 6 close to 
the zero-temperature critical lineH]. This implies 
that the thermal line turns into a line of a first- 
order phase transition while running close to the 
tip of the cusp of the finite-temperature critical 
line for these cases. 

For Nf > 7 the QCDPAX Collaboration re- 
ported that the critical line disappears even in the 
limit of strong coupling. Understanding this phe- 
nomenon by an extension of our line of analysis 



is an interesting problem left for future studies. 

4.4. Continuum limit 

When we increase the temporal lattice size Nt, 
we expect both the thermal line and the cusp of 
the critical line to move toward larger (3. The dis- 
tance between the thermal line and the tip of the 
cusp will diminish, probably as 0{a) « 0{1/Nt). 
To elucidate the properties of chiral transition in 
the continuum limit, one should measure ther- 
modynamic observables across the thermal line 
in the region around the tip of the cusp and ex- 
amine how they vary as a function of Nf. 

An important point to note is that the Nf — 2 
transition is a continuous crossover for finite val- 
ues of Nt- A second-order chiral phase transition, 
as suggested by continuum sigma model analyses 
for this case, would emerge only in the continuum 
limit Nt oo. 
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